Vittrup Man-The life-history of a genetic foreigner in Neolithic Denmark.

The lethally maltreated body of Vittrup Man was deposited in a Danish bog, probably as part of a ritualised sacrifice. It happened between c. 3300 and 3100 cal years BC, i.e., during the period of the local farming-based Funnel Beaker Culture. In terms of skull morphological features, he differs from the majority of the contemporaneous farmers found in Denmark, and associates with hunter-gatherers, who inhabited Scandinavia during the previous millennia. His skeletal remains were selected for transdisciplinary analysis to reveal his life-history in terms of a population historical perspective. We report the combined results of an integrated set of genetic, isotopic, physical anthropological and archaeological analytical approaches. Strontium signature suggests a foreign birthplace that could be in Norway or Sweden. In addition, enamel oxygen isotope values indicate that as a child he lived in a colder climate, i.e., to the north of the regions inhabited by farmers. Genomic data in fact demonstrates that he is closely related to Mesolithic humans known from Norway and Sweden. Moreover, dietary stable isotope analyses on enamel and bone collagen demonstrate a fisher-hunter way of life in his childhood and a diet typical of farmers later on. Such a variable life-history is also reflected by proteomic analysis of hardened organic deposits on his teeth, indicating the consumption of forager food (seal, whale and marine fish) as well as farmer food (sheep/goat). From a dietary isotopic transect of one of his teeth it is shown that his transfer between societies of foragers and farmers took place near to the end of his teenage years.

100 ancient genomes show repeated population turnovers in Neolithic Denmark.

Major migration events in Holocene Eurasia have been characterized genetically at broad regional scales1-4. However, insights into the population dynamics in the contact zones are hampered by a lack of ancient genomic data sampled at high spatiotemporal resolution5-7. Here, to address this, we analysed shotgun-sequenced genomes from 100 skeletons spanning 7,300 years of the Mesolithic period, Neolithic period and Early Bronze Age in Denmark and integrated these with proxies for diet (13C and 15N content), mobility (87Sr/86Sr ratio) and vegetation cover (pollen). We observe that Danish Mesolithic individuals of the Maglemose, Kongemose and Ertebølle cultures form a distinct genetic cluster related to other Western European hunter-gatherers. Despite shifts in material culture they displayed genetic homogeneity from around 10,500 to 5,900 calibrated years before present, when Neolithic farmers with Anatolian-derived ancestry arrived. Although the Neolithic transition was delayed by more than a millennium relative to Central Europe, it was very abrupt and resulted in a population turnover with limited genetic contribution from local hunter-gatherers. The succeeding Neolithic population, associated with the Funnel Beaker culture, persisted for only about 1,000 years before immigrants with eastern Steppe-derived ancestry arrived. This second and equally rapid population replacement gave rise to the Single Grave culture with an ancestry profile more similar to present-day Danes. In our multiproxy dataset, these major demographic events are manifested as parallel shifts in genotype, phenotype, diet and land use.

Population genomics of post-glacial western Eurasia.

Western Eurasia witnessed several large-scale human migrations during the Holocene1-5. Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations.

Medium dependent optical response in ultra-fine plasmonic nanoparticles.

We study the influence of media on the interaction of ultra-fine plasmonic nanoparticles (≤ 8 nm) with radiation. The important role of the surface layer of the nanoparticles, with properties that differ from the ones in the inner part, is established. Using an atomistic representation of the nanoparticle material and its interaction with light, we find a highly inhomogeneous distribution of the electric field inside and around the particles. It is predicted that with an increase in the refractive index of the ambient medium, the extension of the surface layer of atoms increases, something that also is accompanied by an enhanced red shift of the plasmon resonance band compared to large particles in which the influence of this layer and its relative volume is reduced. It is shown that the physical origin for the formation of a surface layer of atoms near the nanoparticle boundary is related to the anisotropy of the local environment of atoms in this layer which changes the conditions for the interaction of neighboring atoms with each other and with the incident radiation. It is shown that a growth of the refractive index of the ambient medium results in an increase in the local field in the dielectric cavity in which a plasmonic nanoparticle is embedded and which is accompanied by a growth of the amplitude of the plasmon resonance. We predict that in the ultra-fine regime the refractive index sensitivity shows a decreasing trend with respect to size which is opposite to that for larger particles. With the applied atomistic model this work demonstrates close relations between field distributions and properties of ultra-fine nanoparticles.

Thermal degradation of optical resonances in plasmonic nanoparticles.

The dependence of plasmon resonance excitations in ultrafine (3-7 nm) gold nanoparticles on heating and melting is investigated. An integrated approach is adopted, where molecular dynamics simulations of the spatial and temporal development of the atoms constituting the nanoparticles generate trajectories out of which system conformations are sampled and extracted for calculations of plasmonic excitation cross sections which then are averaged over the sample configurations for the final result. The calculations of the plasmonic excitations, which take into account the temperature- and size-dependent relaxation of the plasmons, are carried out with a newly developed Extended Discrete Interaction Model (Ex-DIM) and complemented by multilayered Mie theory. The integrated approach clearly demonstrates the conditions for suppression of the plasmons starting at temperatures well below the melting point. We have found a strong inhomogeneous dependence of the atom mobility in the particle crystal lattice increasing from the center to its surface upon the temperature growth. The plasmon resonance suppression is associated with an increase of the mobility and in the amplitude of phonon vibrations of the lattice atoms accompanied by electron-phonon scattering. This leads to an increase in the relaxation constant impeding the plasmon excitation as the major source of the suppression, while the direct contribution from the increase in the lattice constant and its chaotization at melting is found to be minor. Experimental verification of the suppression of surface plasmon resonance is demonstrated for gold nanoparticles on a quartz substrate heated up to the melting temperature and above.

Genomic Steppe ancestry in skeletons from the Neolithic Single Grave Culture in Denmark.

The Gjerrild burial provides the largest and best-preserved assemblage of human skeletal material presently known from the Single Grave Culture (SGC) in Denmark. For generations it has been debated among archaeologists if the appearance of this archaeological complex represents a continuation of the previous Neolithic communities, or was facilitated by incoming migrants. We sampled and analysed five skeletons from the Gjerrild cist, buried over a period of c. 300 years, 2600/2500-2200 cal BCE. Despite poor DNA preservation, we managed to sequence the genome (>1X) of one individual and the partial genomes (0.007X and 0.02X) of another two individuals. Our genetic data document a female (Gjerrild 1) and two males (Gjerrild 5 + 8), harbouring typical Neolithic K2a and HV0 mtDNA haplogroups, but also a rare basal variant of the R1b1 Y-chromosomal haplogroup. Genome-wide analyses demonstrate that these people had a significant Yamnaya-derived (i.e. steppe) ancestry component and a close genetic resemblance to the Corded Ware (and related) groups that were present in large parts of Northern and Central Europe at the time. Assuming that the Gjerrild skeletons are genetically representative of the population of the SGC in broader terms, the transition from the local Neolithic Funnel Beaker Culture (TRB) to SGC is not characterized by demographic continuity. Rather, the emergence of SGC in Denmark was part of the Late Neolithic and Early Bronze Age population expansion that swept across the European continent in the 3rd millennium BCE, resulting in various degrees of genetic replacement and admixture processes with previous Neolithic populations.

Atomistic description of plasmonic generation in alloys and core shell nanoparticles.

Using the extended discrete interaction model we investigate the tunabilty of surface plasmon resonances in alloys and core-shell nanoparticles made from silver and gold in the small (1-15 nm) nanoscale regime where classical models based on the bulk dielectric constant may not apply. We show that the surface plasmon resonance of these alloys and core-shell particles to a large extent follow Vegard's law irrespective of the geometry of the nanoparticle. The evolution of the polarizability with size demonstrates a highly non-linear behaviour of the polarizability with the ratio of the constituents and geometry in alloys and core-shell nanoparticles, with the exception of the longitudinal surface plasmon resonance in nanorods and, partly, nanodisc alloys. We here show that the non-linear behaviour can be explained in terms of the difference in polarizability of the mixing constituents and local effects causing a quenching of the dipoles for geometries with a low aspect ratio. A thorough statistical investigation reveals that there is only a small dependence of the surface plasmon resonance on atomic arrangement and exact distribution in a nanoparticle and that the standard deviation decreases rapidly with the size of the nanoparticles. The physical ground for the random distribution algorithm for alloys in discrete interaction models is explained in detail and verified by a statistical analysis. For nanoparticles below 4 nm a sampling strategy is recommended.

An integrated analysis of Maglemose bone points reframes the Early Mesolithic of Southern Scandinavia.

The extensive peat bogs of Southern Scandinavia have yielded rich Mesolithic archaeological assemblages, with one of the most iconic artefacts being the bone point. Although great in number they remain understudied. Here we present a combined investigation of the typology, protein-based species composition, and absolute chronology of Maglemosian bone points. The majority of the bone points are made from cervids and bovines. However, changes both in species composition and barb morphology can be directly linked to a paucity of finds lasting nearly 600 years in Southern Scandinavia around 10,300 cal BP. We hypothesize that this hiatus was climate-driven and forced hunter-gatherers to abandon the lakes. Furthermore, the marked change in bone points coincides with a change in lithic technology. We, therefore, propose that the Maglemose culture in Southern Scandinavia is fundamentally divided into an Early Complex and a Late Complex.

Exact semi-classical light-matter interaction operator applied to two-photon processes with strong relativistic effects.

X-ray processes involve interactions with high-energy photons. For these short wavelengths, the perturbing field cannot be treated as constant, and there is a need to go beyond the electric-dipole approximation. The exact semi-classical light-matter interaction operator offers several advantages compared to the multipole expansion such as improved stability and ease of implementation. Here, the exact operator is used to model x-ray scattering in metal K pre-edges. This is a relativistic two-photon process where absorption is dominated by electric-dipole forbidden transitions. With the restricted active space state-interaction approach, spectra can be calculated even for the multiconfigurational wavefunctions including second-order perturbation. However, as the operator itself depends on the transition energy, the cost for evaluating integrals for hundreds of thousands unique transitions becomes a bottleneck. Here, this is solved by calculating the integrals in a molecular-orbital basis that only runs over the active space, combined with a grouping scheme where the operator is the same for close-lying transitions. This speeds up the calculations of single-photon processes and is critical for the modeling of two-photon scattering processes. The new scheme is used to model Kα resonant inelastic x-ray scattering of iron-porphyrin complexes with relevance to studies of heme enzymes, for which the total computational time is reduced by several orders of magnitude with an effect on transition intensities of 0.1% or less.

Plasmonic nano-shells: atomistic discrete interaction versus classic electrodynamics models.

Using the extended discrete interaction model and Mie theory, we investigate the tunability of the optical polarizability of small metallic nano-shells. We show that the spectral positions of symmetric and antisymmetric dipolar plasmon resonances vary with the ratio of particle radius to hole radius in a manner similar to one predicted for uniform metallic nano-shells using a semiclassical approach of two coupled harmonic oscillators. We show that, according to the extended discrete interaction model, the dipolar plasmon resonances are also present for nano-shells in the 2-13 nm size region and show the same functional dependence seen for larger nano-shells. Using previously fitted data from experiment, we can predict the size-dependence of the plasma frequency for nano-shells in the 1-15 nm size region. We find that Mie theory, which utilizes the electron mean free path correction for the permittivity, is not able to reproduce the same functional form of the dipolar modes for the nano-shells of the same sizes.

Setting the stage for theoretical x-ray spectra of the H2S molecule with multi-configurational quantum chemical calculations of the energy landscape.

In the H2S molecule, the interplay between different core levels can be investigated in great detail in relation to x-ray spectroscopy, which requires a theory for interpretation. Hence, valence and core excitations into the two antibonding molecular orbitals of the H2S molecule have been calculated within a multi-configurational wave function framework. Scanning along the S-H stretching coordinates, we derive potential energy surfaces and transition dipole moments involving the ground state and core and valence excited states. Both valence excitations and the S1s-1 and S2p-1 core excitations show pairs of dissociative and bound electronic states. These pairs of states are nearly degenerate in H2S at the ground state geometry. The close degeneracy together with conical intersections makes H2S an interesting target for x-ray spectroscopy involving ultra-fast dissociation influenced by non-adiabatic transitions and interference. For future investigations with x-ray absorption spectroscopy (XAS) and resonant inelastic x-ray scattering (RIXS), it is valuable to compare H2S with the water molecule, which exhibits state-selective gating to different vibrational modes [R. C. Couto et al., Nat. Commun. 8, 14165 (2017)] in its well-separated O1s-1 core excited states. The dense manifolds of the S2p-1 core excited states will complicate the analysis of Kα edge RIXS, but dynamical effects could be evaluated through detuning and by comparing with L edge XAS. In L edge RIXS, the dynamical effects will be more pronounced due to the longer lifetime of the S2p-1 core excited states compared to the S1s-1 core excited states.

OpenMolcas: From Source Code to Insight.

In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

Nordic Seaweed and Diabetes Prevention: Exploratory Studies in KK-Ay Mice.

BACKGROUND: The global epidemic of type 2 diabetes (T2D) is a challenging health problem. Lifestyle changes, including nutrition therapy, areimportant for the prevention and management of T2D. Seaweeds contain several bioactive substances with potential health properties and may be a low-cost alternative functional food in the prevention of T2D. OBJECTIVE: The aim of this study was to explore the preventive effects of dried Nordic seaweed species on diabetes in an animal model of T2D. METHOD: Fiftymale KK-Ay mice were randomly assigned to one of four diets: control diet (chow) or diets supplemented with Alaria esculenta (AE), Saccharina latissima (SL), or Palmaria palmata (PP). The effect of the interventions on the progression of T2D was monitored over 10 weeks and evaluated by circulating glucose, glycated hemoglobin (HbA1c), insulin, glucagon, and lipid levels. RESULTS: The SL group had significantly lower bodyweight, lower HbA1c and insulin levels, as well as higher high density lipoprotein (HDL) cholesterol levels after the 10-week intervention than the control group. At the end of the study, the control group had significantly higher HbA1c (p < 0.001) than all of the seaweed groups. CONCLUSION: All seaweed groups improved HbA1C compared to control and Saccharinalatissima seaweed had concomitantly beneficial effects on glycemic control and lipid levels in KK-Ay diabetic mice.

Efficient calculations of a large number of highly excited states for multiconfigurational wavefunctions.

Electronically excited states play important roles in many chemical reactions and spectroscopic techniques. In quantum chemistry, a common technique to solve excited states is the multiroot Davidson algorithm, but it is not designed for processes like X-ray spectroscopy that involves hundreds of highly excited states. We show how the use of a restricted active space wavefunction together with a projection operator to remove low-lying electronic states offers an efficient way to reach single and double-core-hole states. Additionally, several improvements to the stability and efficiency of the configuration interaction (CI) algorithm for a large number of states are suggested. When applied to a series of transition metal complexes the new CI algorithm does not only resolve divergence issues but also leads to typical reduction in computational time by 70%, with the largest savings for small molecules and large active spaces. Together, the projection operator and the improved CI algorithm now make it possible to simulate a wide range of single- and two-photon spectroscopies. © 2019 Wiley Periodicals, Inc.

Fingerprinting Electronic Structure of Heme Iron by Ab Initio Modeling of Metal L-Edge X-ray Absorption Spectra.

The capability of the multiconfigurational restricted active space approach to identify electronic structure from spectral fingerprints is explored by applying it to iron L-edge X-ray absorption spectroscopy (XAS) of three heme systems that represent the limiting descriptions of iron in the Fe-O2 bond, ferrous and ferric [Fe(P)(ImH)2]0/1+ (P = porphine, ImH = imidazole), and FeII(P). The level of agreement between experimental and simulated spectral shapes is calculated using the cosine similarity, which gives a quantitative and unbiased assignment. Further dimensions in fingerprinting are obtained from the L-edge branching ratio, the integrated absorption intensity, and the edge position. The results show how accurate ab initio simulations of metal L-edge XAS can complement calculations of relative energies to identify unknown species in chemical reactions.

Estimating the relative importance of demic and cultural diffusion in the spread of the Neolithic in Scandinavia.

Using a database of early farming sites in Scandinavia, we estimate that the spread rate of the Neolithic was in the range 0.44-0.66 km yr-1 This is substantially slower (by about 50%) than the rate in continental Europe. We interpret this result in the framework of a new mathematical model that includes horizontal cultural transmission (acculturation), vertical cultural transmission (interbreeding) and demic diffusion (reproduction and dispersal of farmers). To parametrize the model, we estimate reproduction rates of early farmers using archaeological data (sum-calibrated probabilities for the dates of early Neolithic Scandinavian sites) and use them in a wave-of-advance model for the first time. Comparing the model with the archaeological data, we find that the percentage of the spread rate due to cultural diffusion is below 50% (except for very extreme parameter values, and even for them it is below 54%). This strongly suggests that the spread of the Neolithic in Scandinavia was driven mainly by demic diffusion. This conclusion, obtained from archaeological data, agrees qualitatively with the implications of ancient genetic data, but the latter are yet too few in Scandinavia to produce any quantitative percentage for the spread rate due to cultural diffusion. We also find that, on average, fewer than eight hunter-gatherers were incorporated in the Neolithic communities by each group of 10 pioneering farmers, via horizontal and/or vertical cultural transmission.

Fluid dynamics in a full-scale flat sheet MBR, an experimental and numerical study.

Fluid dynamics is used for fouling mitigation in membrane bioreactors (MBRs), whereby a proper understanding of the fluid dynamics is of great interest. The influence of fluid dynamics has led to the use of computational fluid dynamics for optimizing MBR systems. In this work, a model has been validated for flat sheet membranes, with use of the Eulerian multiphase method. The model is validated against a comparable setup where the liquid velocities are measured with a laser Doppler anemometer (LDA). Furthermore, the Eulerian multiphase approach is validated against the more numerical direct volume of fluid (VOF) approach with sludge properties for the liquid, resulting in an error between the models of less than 2% for the wall shear stresses. The VOF model further showed that the horizontal components contribute significantly to the total wall shear stresses. The model has been applied to a full-scale setup for studying the effect of deflecting membranes as deflections have been seen in production. Minimizing the deflection of the membrane sheets was crucial to achieve a good operating condition as a deflection of 2 mm in a setup with a gap of 7 mm decreased the wall shear stresses with as much as 40% on average on the specific membrane surface.

Molecular and Electronic Structure of Re2Br4(PMe3)4.

The dinuclear rhenium(II) complex Re2Br4(PMe3)4 was prepared from the reduction of [Re2Br8](2-) with (n-Bu4N)BH4 in the presence of PMe3 in propanol. The complex was characterized by single-crystal X-ray diffraction (SCXRD) and UV-visible spectroscopy. It crystallizes in the monoclinic C2/c space group and is isostructural with its molybdenum and technetium analogues. The Re-Re distance (2.2521(3) Å) is slightly longer than the one in Re2Cl4(PMe3)4 (2.247(1) Å). The molecular and electronic structure of Re2X4(PMe3)4 (X = Cl, Br) were studied by multiconfigurational quantum chemical methods. The computed ground-state geometry is in excellent agreement with the experimental structure determined by SCXRD. The calculated total bond order (2.75) is consistent with the presence of an electron-rich triple bond and is similar to the one found for Re2Cl4(PMe3)4. The electronic absorption spectrum of Re2Br4(PMe3)4 was recorded in benzene and shows a series of low-intensity bands in the range 10 000-26 000 cm(-1). The absorption bands were assigned based on calculations of the excitation energies with the multireference wave functions followed by second-order perturbation theory using the CASSCF/CASPT2 method. Calculations predict that the lowest energy band corresponds to the δ* → σ* transition, while the next higher energy bands were attributed to the δ* → π*, δ → σ*, and δ → π* transitions.

Molecular Orbital Simulations of Metal 1s2p Resonant Inelastic X-ray Scattering.

For first-row transition metals, high-resolution 3d electronic structure information can be obtained using resonant inelastic X-ray scattering (RIXS). In the hard X-ray region, a K pre-edge (1s→3d) excitation can be followed by monitoring the dipole-allowed Kα (2p→1s) or Kß (3p→1s) emission, processes labeled 1s2p or 1s3p RIXS. Here the restricted active space (RAS) approach, which is a molecular orbital method, is used for the first time to study hard X-ray RIXS processes. This is achieved by including the two sets of core orbitals in different partitions of the active space. Transition intensities are calculated using both first- and second-order expansions of the wave vector, including, but not limited to, electric dipoles and quadrupoles. The accuracy of the approach is tested for 1s2p RIXS of iron hexacyanides [Fe(CN)6](n-) in ferrous and ferric oxidation states. RAS simulations accurately describe the multiplet structures and the role of 2p and 3d spin-orbit coupling on energies and selection rules. Compared to experiment, relative energies of the two [Fe(CN)6](3-) resonances deviate by 0.2 eV in both incident energy and energy transfer directions, and multiplet splittings in [Fe(CN)6](4-) are reproduced within 0.1 eV. These values are similar to what can be expected for valence excitations. The development opens the modeling of hard X-ray scattering processes for both solution catalysts and enzymatic systems.